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Monolayers of transition metal dichalcogenides, namely, 
molybdenum and tungsten disulfides and diselenides 
demonstrate unusual optical properties related to the 
spin-valley locking effect. Particularly, excitation of 
monolayers by circularly polarized light selectively cre¬ 
ates electron-hole pairs or excitons in non-equivalent 
valleys in momentum space, depending on the light Fe¬ 
licity. This allows studying the inter-valley dynamics of 
charge carriers and Coulomb complexes by means of op¬ 
tical spectroscopy. Here we present a concise review of 
the neutral exciton fine structure and its spin and val¬ 
ley dynamics in monolayers of transition metal dichalco¬ 
genides. It is demonstrated that the long-range exchange 
interaction between an electron and a hole in the exci¬ 
ton is an efficient mechanism for rapid mixing between 
bright excitons made of electron-hole pairs in different 
valleys. We discuss the physical origin of the long-range 
exchange interaction and outline its derivation in both 
the electrodynamical and k • p approaches. We further 
present a model of bright exciton spin dynamics driven 
by an interplay between the long-range exchange inter¬ 
action and scattering. Finally, we discuss the application 
of the model to describe recent experimental data ob¬ 
tained by time-resolved photoluminescence and Kerr ro¬ 
tation techniques. 
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Kerr rotation dynamics at T = 4 K for a cr"'" and a~ pump 
beam in WSe 2 . Inset: Schematics of the optical selection rules 
of the excitons photogenerated from charge carriers in K± val¬ 
leys and their coupling induced by the long-range exchange in¬ 
teraction. From [C. R. Zhu, K. Zhang, M. Glazov, B. Urbaszek, 
T. Amand, Z. W. Ji, B. L. Liu, and X. Marie, Exciton valley 
dynamics probed by Kerr rotation in WSe2 monolayers, Phys. 
Rev. B 90, 161302(R) (2014)]. 
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1 Introduction Transition metal dichalcogenide mono- 
layers have remarkable electronic, optical, optoelectronic 
and chemical properties [1-3]. They can serve as a build¬ 


ing blocks for novel type of heterostructures, known as 
van der Waals heterostructures, where layers of different 
materials can be isolated and then assembled in a tai- 
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lored sequence [4]. Monolayers of MX 2 where M=Mo,W 
stands for the transition metal and X=S,Se,Te stands for 
the chalcogen have a hexagonal lattice like graphene. In 
contrast to graphene, where conduction and valence bands 
merge at the edges of the Brillouin zone [5], they demon¬ 
strate direct bandgaps ^ 2 eV realized at the edges, points 
K±, of the Brillouin zone [3]. 

Strong spin-orbit interaction in MX 2 monolayers leads 
to prominent effects termed as spin-valley locking [6]. 
First, valence and conduction band states are spin split 
already at the K± points [7-11]. The splitting is of op¬ 
posite signs in the and K_ valleys and acts as a 
valley-dependent static magnetic field applied perpendic¬ 
ular to the monolayer plane. Second, optical transition 
rules for interband transitions in the valleys K± are chiral: 
cr+ light induces transitions in one valley of the energy 
spectrum, while a~ light induces transitions in the 
other one, K_ [12-14]. Hence, excitation of MX 2 sam¬ 
ples by circularly polarized light induces both spin and 
valley polarization of charge carriers. These phenomena 
make transition metal dichalcogenides ideal systems for 
studying spin and valley dynamics in both experiment and 
theory, see for example Refs. [6,12-17]. 

Substantial optical orientation in transition metal 
dichalcogenides such as M 0 S 2 and WSe 2 related to selec¬ 
tive excitation of the valleys by circularly polarized light 
has been indeed revealed in experiments, e.g. [16,18-24]. 
Moreover, single particle spin states and valley coherence 
were expected to be extremely robust for several inter¬ 
linked reasons: First, it is necessary to transfer a huge 
wavevector, on the order of inverse lattice constant, for the 
charge carrier to change valley. Second, the carrier needs 
also to flip its spin, which is energetically unfavourable 
due to the large spin-orbit splittings [25-29]. Despite the 
valley index stability expected in this single particle pic¬ 
ture, recent optical spectroscopy experiments demonstrate 
surprisingly rapid intervalley transfer on a picosecond time 
scale [16,22-24]. 

In reality the optical properties of MX 2 monolayers are 
governed not by free carriers but by excitons, electron-hole 
pairs strongly bound by Coulomb attraction. In this fea¬ 
ture article we will show in detail why and how excitons 
govern the polarization dynamics observed in optical spec¬ 
troscopy. This allows understanding how the optically gen¬ 
erated spin/valley polarization can decay on a picosecond 
timescale, in good agreement with the experimental find¬ 
ings [30]. Indeed, weak screening of the Coulomb interac¬ 
tion and relatively large effective masses of electrons and 
holes result in strong, tightly-bound, excitons in transition 
metal dichalcogenides monolayers. The exciton binding 
energies in these materials ^ 0.5 eV [31-35] exceed by 
far those in conventional two-dimensional systems based 
GaAs-like semiconductors where Efj < 10“^ eV [36-38]. 
Experiments and theoretical modeling reveal also photo¬ 
luminescence lifetimes in the picosecond range indicat¬ 
ing high exciton oscillator strengths of excitons in MX 2 


o o o 

I I • ^ I I 

6 ..., .,.- 6 -,, ...- 6 -,, jb 

""" 0 ""' "" 0 '^‘ "" 0 "'’ 

Figure 1 Schematic illustration of atoms in hexagonal lat¬ 
tice of MX 2 . Yellow circles denote positions of transition 
metal M, blue circles show the projections on the mono- 
layer plane of the chalcogens X; 0 denotes the origin of 
point transformations, Ri, R 2 and Rs are the position vec¬ 
tors of metal atoms. 


monolayers [20,39-41]. Despite strong suppression of in¬ 
dividual carrier spin flips, the exchange interaction in the 
electron-hole pair forming an exciton [30, 42, 43] turns out 
to be strong enough to enable simultaneous spin flip of 
an electron and a hole [44]. The exciton spin dynamics 
in monolayers of transition metal dichalcogenides is, thus, 
similar to that of excitons in conventional two-dimensional 
semiconductor structures [38]. 

Here we provide a brief overview of exciton fine struc¬ 
ture and spin dynamics in monolayers MX 2 . The article 
is organized as follows: We start in Sec. 2 with the ba¬ 
sic introduction to the band structure of transition metal 
dichalcogenides monolayers and to excitonic effects in 
these materials. A brief overview of exciton complexes 
investigated in theory and experiment in MX 2 monolayers 
is given in Sec. 3. Further we discuss the fine structure of 
the bright excitons and review theoretical approaches to 
calculate the exchange interaction between an electron and 
a hole in Sec. 4. Next, in Sec. 5 we address theoretically 
exciton spin dynamics governed by the exchange interac¬ 
tion in MX 2 monolayers and discuss recent experimental 
findings on excitonic spin dynamics in Sec. 6. The outlook 
and conclusions are presented in Sec. 7 

2 Band structure The band structure of transition 
metal dichalcogenide monolayers is studied theoretically 
by various methods, see Ref. [45] for a review, including 
both first principles approaches, e.g. [10,26,31,46,47], and 
empirical ones. The latter involve empirical tight-binding 
models [9, 48-50] and the k • p perturbation theory [7, 8, 
11,25]. For our purposes it is enough to use the simplest 
possible, two-band (or four-band, if spin is included) k • p 
method [7]. To establish the form of effective k p Hamilto¬ 
nian we make use of symmetry arguments and the method 
of invariants. 

Figure 1 represents the positions of atoms in the hexag¬ 
onal lattice of MX 2 projected onto the monolayer plane. 
The symmetry of the monolayer is described by the point 
group Dsh which includes the horizontal refiection plane 
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containing metals, and three-fold rotation axis C 3 in¬ 
tersecting the horizontal plane in the center of hexagon. 
The complete list of Ds/^ elements includes also the mirror- 
rotation axis S's, three two-fold rotation axes C 2 lying in 
the monolayer plane and mirror reflection planes cFy which 
contain the C 2 axes. 

The Brillouin zone of MX 2 monolayers is also hexag¬ 
onal and the direct band gaps are realized at the points K± 
corresponding to its edges, see inset in Fig. 2. These two 
valleys are coupled by the time-reversal symmetry. The 
wave vector group at the K± points is The symme¬ 
try of conduction and valence band wavefunctions can be 
established taking into account that the conduction band 
is formed, neglecting spin and spin-orbit coupling, mostly 
from the orbitals of metal, while the valence band is 
formed mostly from the orbitals [7,47]. Here x 

and y are the Cartesian axes in the monolayer plane and 2 ; 
is the sample normal. The state symmetry at the F point is 
completely determined by that of the atomic orbitals. By 
contrast, the symmetry behavior of the states at the edges 
K± of the Brillouin zone is determined not only by the 
symmetry of the particular orbitals but also by the rela¬ 
tion between the phase K± • Rj of the Bloch function and 

the phase (^CsK±^ • Rj, where Rj is the orbital position 

and Cs is the three-fold rotation operator. As a result [ 8 ], 
the valence band orbital Bloch functions being oc 
in the and K- valleys, respectively, transforms ac¬ 
cording to the invariant representation of the Csh point 
group. A' in notations of Ref. [ 8 ] or Fi in the notations 
of Ref. [51]. The conduction band states transform in the 
valleys and K- according to E[ (F 2 ) like x + iy, and 
E 2 (Is), like X — iy, respectively. This symmetry analysis 
demonstrates that selection rules for the interband optical 
transition at the normal incidence of radiation are chiral: 
The transitions in the valley are active in the cr+ po¬ 
larization of light, while the transitions at K- are excited 
by the (j~ polarized light [12-14]. 

Inclusion of the electron spin and spin-orbit coupling 
does not change these chiral selection rules. Indeed, tak¬ 
ing into account the electron spin, the valence band states 
transform according to the spinor representations Iy, spin- 
up, and Ig, spin-down, in each valley, see Fig. 2. The con¬ 
duction band states in the Kj^ valley transform according 
to Fii, spin-up, and Ig, spin-down, while the K_ valley 
conduction band states transform according to Fio, spin- 
up, and Fi 2 , spin-down, respectively. These spinor repre¬ 
sentations of the Cs/j, point group are non-degenerate [51] 
which means that, in the absence of external forces, both 
the conduction and valence bands already are spin split 
at the K± points. The time-reversal symmetry ensures 
the splittings to be of opposite signs in the Kj^ and K_ 
valleys. Therefore, as shown in Fig. 2, the order of spin 
states at the and K- edges of the Brillouin zone is 
reversed. The spin splittings in each valley are sizeable: 
For the valence band states their order of magnitude is 



Figure 2 Sketch of the band structure of transition metal 
dichalcogenide monolayers. The bands are labelled by the 
corresponding irreducible spinor representations with ar¬ 
rows in parentheses indicating the electron spin orienta¬ 
tion. Solid and dashed arrows show the transitions active 
under the normal incidence in the cr+ and cf ~ polarizations, 
respectively. Inset shows the Brillouin zone and the band 
dispersion. The order of conduction band states is shown 
for M 0 S 2 in accordance with Ref. [9]. 


~ 100 meV increasing with the atomic number of con¬ 
stitutive elements, while for the conduction band the split¬ 
tings are in the meV to tens of meV range [-1 ]. The 
dipole interaction with electromagnetic held conserves the 
spin, therefore, there are only two allowed interband opti¬ 
cal transitions under the normal incidence at each valley, as 
illustrated in Fig. 2 by vertical arrows. The selection rules 
are determined by the orbital part of the Bloch functions. 
Hence, the transitions in the K± valleys are activated by 
-polarized photons, respectively. 

In what follows we focus on the optical transitions in¬ 
volving the bottom conduction and topmost valence bands, 
i.e., the states TV and Ai in the valley and the states Eg 
and Fi 2 in the K- valley. Therefore, it is enough to resort 
to the two-band model where the effective k • p Hamiltoni¬ 
ans in the vicinity of K± points can be presented as 2 x 2 
matrices, respectively. 


%+ 


( Eg j3k-\ 
V73fc+ 0 


H- = 


( Eg -73fc+\ 

\-l 3 k- 0 j ' 

( 1 ) 


Here the energy is referred to the valence band top. Eg is 
the band gap, k± = kx E iky, k = {kx^ky) is a two- 
dimensional wave vector of the electron reckoned from the 
corresponding Brillouin zone edge, and the real constant 
73 is related with the interband matrix element of the elec¬ 
tron momentum. Note that the choice of signs in the off- 
diagonal terms in'H± as well as the fact that 73 parameter 
is real is related with a deflnite convention about the phases 
of the Bloch functions and action of the time-reversal op¬ 
erator. In the two-band approximation neglecting contri¬ 
butions from distant bands and from bare electron disper- 
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sion [52] the conduction and valence bands are symmetric 
with the effective masses rric = —rriy = m*, where 


m* 


( 2 ) 


In MX 2 monolayers the effective mass m* is about 0.5mo, 
where mo is the free electron mass. More elaborate k • p 
models are presented in Refs. [7, 8 ,11, 52, 53] and include 
remote bands into the consideration leading to the electron- 
hole asymmetry and a better agreement with atomistic cal¬ 
culations. For instance, Ref. [53] reports me = 0.56mo, 
my = —0.59mo for MoSe 2 and my = 0.28mo, my = 
—0.36mo for WSe 2 . It shows that the above two-band ap¬ 
proximation can be used if a moderate accuracy for effec¬ 
tive masses is needed. Overall, existing k • p models [53] 
provide adequate results for the dispersion of electrons and 
holes in the vicinity of K± points, but do not reproduce 
the observed Zeeman splitting for excitons and trions [52]. 
However, the fine structure of the radiative doublet of the 
excitonic states, which is the focus of this paper, is not af¬ 
fected by the details of electron and hole bands. Hence, 
hereafter we use the simplest possible band structure mod¬ 
els that will be particularly useful to describe the valley 
polarization dynamics. 


3 Exciton states in MX 2 monolayers Relatively 
heavy masses of charge carriers and weak screening in 
two-dimensional systems yield strong excitonic effects, 
see Ref. [54] for recent review. It is worth mentioning, that 
strong Coulomb effects in MX 2 monolayers make it pos¬ 
sible to observe, besides neutral excitons, other Coulomb 
correlated complexes such as trions (also known as charged 
excitons) [32] and biexcitons [63]. The spin and valley dy¬ 
namics of these complexes deserves further experimental 
and theoretical study and is not addressed here. 

In the following we assume that the main contribution 
to the exciton wavefunctions is provided by the conduc¬ 
tion and valence band states in the vicinity of K± points 
of the /c-space and denote the conduction (valence) band 
branches by the spin index Sy (sh) and the valley index Ty 
(Th). Then the exciton state in a transition metal dichalco¬ 
genide monolayer can be presented within the k • p formal¬ 
ism as [64] 

^n]SeTe]Shrh {ke,kh)\SeTeke-,ShThkh), (3) 

^ei^h 

where n = Is, 2 p,... denotes the states of electron-hole 
relative motion, ky (kh) is the electron (hole) wavevec- 
tor, and \syTyky] ShThkh) denotes the excited state of the 
MX 2 monolayer with the occupied conduction band state 
IsyTyky) 311(1 the empty valence band state IC\shThkh), 
with 1C being the time-reversal operator. It is convenient 
to present this two-particle excited state in the coordinate 
representation as a product [64,65] 

ISeTefee; ShThkh) = (rh), (4) 



Energy (eV) 

Figure 3 The derivative of the reflectance contrast spec¬ 
trum, d{AR/R)/duj, of the WS 2 monolayer. The exciton 
ground state and the higher excited states are labeled by 
their respective quantum numbers (schematically shown at 
bottom-right corner). The spectral region around the Is 
transition (AX) and the trion peak (AXt) of the A exci¬ 
ton is scaled by a factor of 0.03 for clarity. Inset shows the 
as-measured reflectance contrast ARjR, allowing for the 
identification of the A, B and C transitions. From Ref. [33]. 


where (^e) and are the conduction and 

valence band Bloch functions, the latter taken in the hole 
representation. As functions of ky and kh, the coefficients 
Cn-,SeTe;shrh{^ej kh) are the Fourier transforms of exciton 
real-space envelope functions Cn-seTe^shThi'^e, rh). 

Neglecting the band mixing the exciton states can be 
labelled by the center of mass wavevector K = (AT^,, Ky), 
the envelope index n as well as by the spin and valley in¬ 
dices of electrons and holes. The substantial spin splitting 
of the valence band allows us to consider independently 
A- and B-excitons formed of holes, respectively, in the 
top/bottom spin valence subbands, see Fig. 2. In what fol¬ 
lows we focus on the A-exciton states. 

Estimates in the framework of Wannier-Mott model 
of two-dimensional exciton give binding energies A 5 = 
2 \ie^ ! (x/i)^ r\j 500 meV for the reduced mass 11 = 
771 * /2 = 0.25mo and the background dielectric constant 
X = 5. The latter follows from the assumption that the 
two-dimensional layer is deposited on a semi-infinite sub¬ 
strate with the static dielectric constant 10. An accurate 
evaluation of the exciton binding energy is a very complex 
problem, mainly due to non-trivial screening effects and 
the necessity to account for fine details of the band struc¬ 
ture. Various theoretical approaches have already been 
developed see, e.g., [47,58-62], which report values of A 5 
in the range of ^ 100 meV to ^ 1 eV depending on the 
material, dielectric surrounding and level of approxima- 
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tions used. Recent works based on reflectivity measure¬ 
ments [33] and two-photon absorption [34] indeed report 
the A-exciton binding energies in the range of above es¬ 
timation for transition metal dichalcogenides monolayers. 
Figure 3 shows the derivative of the reflectance contrast 
spectrum, d{AR/R)/duj of a WS 2 monolayer measured 
in Ref. [33]. Here AR is the difference of the reflection 
coefficients of the sample and substrate, R is the sub¬ 
strate reflection coefficient, uj is the frequency of the light. 
The data clearly demonstrate s-shell excitonic series up 
to the principle quantum number 5. The analysis carried 
out in Ref. [33] reveals strong deviations from the two- 
dimensional hydrogen-like model which is attributed to 
the complex screening of the Coulomb interaction in the 
MX 2 monolayers. The binding energy of exciton in WS 2 
was found to be about 320 meV in close agreement with 
theoretical estimates [33]. The higher value of exciton 
binding energy Ei, ^ 700 meV for the same material was 
extracted in Ref. [56]. A value of Efj ^ 600 meV, was 
reported in [34] for WSe 2 from identiflcation of 25 and 
2p excitonic states by means of single- and two-photon 
spectroscopy and DFT-GW band gap calculation, though a 
smaller value of ^ 370 meV have been also reported [55]. 
A similar value of the binding energy, Ei^ > 570 meV, was 
reported for excitons in M 0 S 2 from the photocurrent spec¬ 
troscopy [57]. The spread of reported binding energies can 
be related with use of different experimental and sample 
preparation methods as well as with unprecise knowledge 
of the single-particle band gap in MX 2 systems. Impor¬ 
tantly, all experimental and theoretical results clearly show 
that the optical properties are governed by robust excitons. 

Now we proceed to the flne structure of excitonic states 
which is the main goal of the paper. For a given rela¬ 
tive motion quantum number n {n = I 5 in what follows) 
and center of mass wavevector K we have eight A-exciton 
states in total, four of them are intra-valley (or direct) ex¬ 
citons with the occupied state in the conduction band and 
the empty state in the valence band being in the same val¬ 
ley (note that in this case Tq = —Th since Th refers to the 
hole valley index) and the others four are inter-valley (indi¬ 
rect) excitons. The latter can manifest themselves in opti¬ 
cal transitions if a third particle, a phonon or an impurity, is 
involved in the transition. Moreover, it follows from the se¬ 
lection rules illustrated in Fig. 2 that in the dipole approxi¬ 
mation only two of four intra-valley excitons are coupled to 
the light under normal incidence: those are the states with 
5e = —Sh (again, spins of empty state in the valence band 
and in the conduction band are the same). The remaining 
two states are dark at the normal incidence. Note, that the 
A-excitons with Se = Sh (spin-forbidden) are, due to the 
spin-orbit mixing, active only in the z-polarization [30]. 
In spite of being dark, these states can play an important 
role in the photoluminescence formation [66] . Other dark 
states, such as 2p excitons are visible in two-photon ab¬ 
sorption processes [34,56]. These effects are beyond the 
scope of the present paper. 



Figure 4 Optical selection rules of the two bright A- 
exciton states with the small center of mass wavevector K 
and their long-range Coulomb exchange coupling, Eq. (5). 



Figure 5 Schematic illustration of exciton dispersion 
showing splitting of the bright doublet into linearly po¬ 
larized eigenstates with the dipole moment oriented par¬ 
allel (L) and perpendicular (T) to the in-plane wavevector. 
K = 0 point corresponds to center of exciton Brillouin 
zone. Red solid lines show photon dispersion (light cone). 
Dispersion is shown not to scale. Inset shows orientation of 
effective held vs. wavevector K direction in Eq. (5), 
the splitting is not shown. 


4 Fine structure of bright excitons Bright A- 
excitonic states in the and K- valleys transform 
according to the 7^2 and Is representations of Csh point 
symmetry group and are excited by cr+ and cf~ photons, 
respectively, see Fig. 4 for illustration. This pair of states 
can be also related with the two dimensional representa¬ 
tion Eq of the Ds/i point group relevant for the overall 
symmetry of the MX 2 monolayer. The symmetry argu¬ 
ments show that for excitons propagating in the structure 
plane with the center of mass wave vector K = {Kx, Ky) 
the eigenstates are linear combinations of the circularly 
polarized states corresponding to the microscopic dipole 
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moment of exciton oscillating parallel or perpendicular to 
its wave vector [38, 64,67, 68]. The effective Hamiltonian 
acting in the basis of / 2 and Is states can be written in the 
framework of invariants method as a 2 x 2 matrix: 

= ^ ■ cr) . (5) 

Here a = a{K) is the coupling parameter to be deter¬ 
mined from the microscopic model. The second line in 
Eq. (5) gives the effective Hamiltonian in the pseudospin 
representation, with a = {(Jx^cFy) being the vector com¬ 
posed of Pauli matrices and is the effective pseudospin 

precession frequency also termed as effective “magnetic- 
like” field acting on an exciton pseudospin 1/2. It has 
non-zero components hf2x = 2aK^ cos 21 ^ and HQy = 
2aK^ sin 2^, where ^ is the angle between K and the in¬ 
plane axis X. The effective field is described by second an¬ 
gular harmonics of the exciton wave vector K, this is be¬ 
cause a transfer of angular momentum 2 is needed to flip 
exciton circular polarization. The dependence of the spin 
precession frequency I7 k on the direction of the exciton 
center of mass wavevector K is schematically presented in 
the inset to Fig. 5. The splitting between the eigenmodes is 
HQk = 2aK ‘^; it corresponds to the exciton longitudinal 
transverse splitting and is depicted in Fig. 5. 

The coupling parameter a{K) in Eq. (5) is determined, 
in fact, by the long-range exchange interaction between an 
electron and a hole [64,65,69]. It can be found either in the 
framework oik p perturbation theory or in the electrody¬ 
namical approach. In the latter case the interaction of the 
mechanical exciton, i.e., the electron-hole pair bound by 
the direct Coulomb interaction [70], with the macroscopic 
longitudinal electric field should be taken into account. We 
start with the electrodynamical treatment and then demon¬ 
strate its equivalence to the /c • p calculation. 

In the first approach the exciton frequencies can be 
found from the poles of the refiection coefficient of the 
two-dimensional structure [70]. Let us derive the refiec¬ 
tion coefficient of the MX 2 monolayer at the frequency 
range close to the lowest in energy excitonic transitions 
(A-excitons). For simplicity we consider the monolayer sit¬ 
uated in {xy) plane surrounded by dielectric media with 
high-frequency a dielectric constant the contrast of 
background dielectric constants between the material and 
surroundings is disregarded. The allowance for the dielec¬ 
tric constant as well as for the substrate is straightforward 
following Ref. [64], see below. The geometry is illustrated 
in Fig. 6. 

Maxwell equations for the electric field vector E can 
be recast as [64] 

rotrotF^=^ —^ + 47rP(z)], (6) 


a{Kx - iKy)^ 


where uj is the frequency of light and P{z) is the exci¬ 
tonic contribution to the dielectric polarization. We solve 
Maxwell equations (6) the for electromagnetic field tak¬ 
ing into account the excitonic contribution to the dielec¬ 
tric polarization. The latter is calculated within the first- 
order perturbation theory using dipole approximation for 
exciton-light coupling and assuming Wannier-Mott exciton 
as [30,64]: 


5{z)\m?E{z) 

ujo — UJ + ir h^ujQ 


(7) 


Here uj is the incident radiation frequency, ujq is the exci¬ 
ton resonance frequency determined by the band gap and 
its binding energy, F is its nonradiative damping, ip{p) is 
the relative electron-hole motion wavefunction in the exci¬ 
ton, and 73 is defined in Eq. (1). The short-range part of 
the exchange interaction contribution can also be included 
in ujQ. Factor 5{z) ensures that the dipole moment is in¬ 
duced only in the monolayer of transition metal dichalco¬ 
genide whose width is negligible compared with the light 
wavelength. At a normal incidence of radiation the ampli¬ 
tude reflection coefficient of a monolayer has a standard 
form [64,71] 


r(w) = 


i^o 

u)o-u) - i(ro + r) ’ 


( 8 ) 


where 


ro = 


27rge^7| 

H^kijUJq 


l^(0)P, 


(9) 


is the radiative decay rate of an exciton in the monolayer, 
Q. = v^cj/c is the wavevector of radiation. Rough es¬ 
timate of radiative lifetime for the Wannier-Mott exciton 
parameters in MX 2 gives l/(2/o) ^ 1 ps in agreement 
with recent measurements [20, 39, 40]. Note, that Iq in 
transition metal dichalcogenides exceeds by far the radia¬ 
tive decay rate for excitons in conventional semiconductor 
quantum wells [36-38,72] The parameters of the pole in 


Q 



X 

Figure 6 Schematic illustration the system geometry with 
s (TE mode) and p (TM mode) polarized incident light. 
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the reflectivity, Eq. (8), describe the eigenenergy and decay 
rate of the exciton with allowance for the light-matter inter¬ 
action. We note that in the approximation of strictly two- 
dimensional layer the mechanical exciton eigenfrequency 
is not renormalized. In this approximation the transmission 
coefficient amplitude is equal to t{uj) = 1 -l-r(cu). In agree¬ 
ment with overall symmetry of the MX 2 monolayer, 
the reflection coefficient at a normal incidence is polariza¬ 
tion independent. 

In order to And the flne structure of the moving exciton 
energy spectrum we consider the oblique incidence of ra¬ 
diation in the (x^)-plane. Fig. 6. The solution of Maxwell 
equations (6) yields two eigenmodes of electromagnetic 
held: s (or TE) polarized wave with ^ and perpen¬ 
dicular to the light incidence plane, and p (TM) polarized 
wave with E in the incidence plane, see Fig. 6. The reflec¬ 
tion coefficients in a given polarization i = s or p have, 
like Eq. (8), the pole contributions with the modifled pa¬ 
rameters 


niui) = 


_ i-Tpi _ 

uioi-ui - i{roi + r) ’ 


( 10 ) 


where 


ro, = ^ro, rop = ^ro, 
Qz q 


( 11 ) 


Qz = {q‘^ — qf\)^^‘^ is the 2 ) component of the light wavevec- 
tor, is its in-plane component, and cuoz = ^o(^||) = 
cuo + /(2M) is the mechanical exciton frequency, M is 
the exciton effective mass. 

The light-matter interaction results in the radiative de¬ 
cay of the excitons with the wavevectors inside the light 
cone, ^ see Fig. 5. The decay rates are differ¬ 

ent for T- and L-polarized excitons. For the excitons out¬ 
side the light cone, qz becomes imaginary and correspond¬ 
ing exciton-induced electromagnetic held decays exponen¬ 
tially with the distance from the monolayer. Therefore, ex¬ 
citon interaction with the held results in renormalization 
of its frequency rather than its decay rate [64, 68, 73]. For¬ 
mally, it corresponds to imaginary F^i in Eqs. (10), (11). 
Making analytical continuation of Eqs. (10) to the complex 
qz, and introducing the notation K = for the center of 
mass wavevector of an exciton, we obtain from the poles 
of reflection coefficients the splitting between the longitu¬ 
dinal (P II K) and transverse (P T K) exciton states: 


AE = hFc 






hFo — , 


( 12 ) 


where the approximate equation holds for iT ^ g and one 
can replace uj by ccq in the deflnition of q since exciton 
dispersion is almost flat as compared with the dispersion 
of light. Equation (12) for excitons outside the light cone is 
in agreement with the phenomenological Hamiltonian (5). 
The comparison of the splittings enables us to determine 
the constant a in Eq. (5) as 


a{K) = 


hFo 

2Kq' 


(13) 


The long-range exchange splitting of excitonic states scales 
approximately (x K for K q, i.e. for the states outside 
the light cone. We stress that the states within the light cone 
undergo renormalization of damping, i.e. imaginary part 
of the energy, rather than of the renormalization of its real 
part. 

Equation (13) can be easily generalized for the realistic 
case of the system “vacuum - monolayer of MX 2 - sub¬ 
strate” under assumption of the same background refrac¬ 
tive index n of the monolayer and the substrate. The reflec¬ 
tion at the boundary “vacuum - monolayer” yields the re¬ 
placement in denominators of Eq. (10) Fq^i FQ^i{l-\-ri), 
where background reflectivities in s and p polarizations 
are given by the Fresnel formula with the refractive index 
n [74]. In this case in accordance with Ref. [23] 


a{K) = 


cHFq n + 1 
2KuJo + 1 


(14) 


Thus, the long-range exchange interaction constant a is 
determined by the radiative decay rate of the exciton and 
the geometry of the system. Therefore, is expected to be 
enhanced in transition metal dichalcogenides monolayers 
compared with conventional semiconductor structures due 
to high oscillator strengths of excitons in MX 2 . 

The electrodynamical treatment is general and is inde¬ 
pendent of the particular model of the band structure and 
of exciton wavefunction: The key parameter in Eqs. (13), 
(14) is the radiative decay rate of exciton, Fq. It can be con¬ 
sidered as a phenomenological parameter of the theory to 
be found from independent experiments. However, it is in¬ 
structive to derive Eqs. (5) and (13) making use of k p per¬ 
turbation theory as well. It allows us to illustrate how the 
Coulomb interaction enables the electron-hole pair transfer 
between the and K- valleys. As a first step, follow¬ 
ing Ref. [30] we consider the exchange interaction between 
two electrons ipm and occupying F 12 band in the K- 
valley and Fj band in the valley with the wave vectors 
ki and k 2 reckoned from the edges of Brillouin zone, see 
red filled circles in Fig. 7(a). Making use of the Hamiltoni¬ 
ans (1) and taking into account k • p interaction in the first 
order we present their wavefunctions in the form 



^9 J 

^n{r2) = 

Upr2) - ^^^U+{r2) 


(15a) 

(15b) 


Here Up{r) are the Bloch amplitudes at the points, 
T = ± enumerates valleys, subscript F denotes the rep¬ 
resentation of Csh point group corresponding to the par¬ 
ticular Bloch function. As the final states of the electron 
pair we consider the conduction band Fn state in the 
valley, 'iprn ', and the valence band Fg state in the K- val¬ 
ley, 2pri' characterized by the wavevectors k[ and 
spectively. Similarly to Eqs. (15) their wavefunctions read. 
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Figure 7 Sketch of exciton intervalley transfer: (a) shows the initial state, (b) shows the final state. Filled circles denote 
electrons, empty circle denotes hole. Dash-dotted arrow shows the Coulomb interaction, green dotted lines show the k • p 
mixing. 


Fig. 7(b), 


c^iKn)- 


73^1.+ rr- 


Er, 




V’n'(r- 2 ) = \u+{r2) - 


(16a) 

(16b) 


According to the general theory [65] the exchange matrix 
element of the Coulomb interaction V (ri—r 2 ) can be writ¬ 
ten as 


{m'n'\V{ri -r 2 )\mn) = -Vk[-kjki+k 2 ,k[+k'^ 

X (17) 

Here Vq is the two-dimensional Fourier transform of the 
bare Coulomb potential. Equation (17) clearly shows that 
the combined action of the k • p mixing and the Coulomb 
interaction enables intervalley transfer of the electron-hole 
excitation, see Fig. 7. As demonstrated in Fig. 7 each elec¬ 
tron remains in the same valley, but due to k ' p mixing 
changes the band. Therefore, the long-range exchange in¬ 
teraction transfers the electron-hole pair between the K_ 
and valleys of the energy spectrum. 

As a second step of the derivation, in order to repro¬ 
duce Eq. (5), we make standard transformations from the 
electron-electron representation to the electron-hole rep¬ 
resentation in Eq. (17), see Eqs. (3), (4) and Ref. [65] 
and perform averaging over the relative motion wavefunc- 
tion. Finally, inclusion of the high-frequency screening, see 
Refs. [68, 75, 76] for details, yields off-diagonal element 
(F'2|7 ^x(K)|/ 3) in Eq. (5). The LT splitting in Eq. (12) in 
agreement with Eq. (9) reads 


with pcv being interband momentum matrix element. We 
stress that the Coulomb interaction is long-range, it does 
not provide intervalley transfer of individual electrons, 
however, the exchange process involves one electron from 
the and another one from the valley. Equations (5) 
and (12) describe the long-range (also termed annihilation 
or resonant) exchange interaction between an electron and 
a hole. It was previously derived for excitons in semicon¬ 
ductor quantum wells both on the basis of electrodynamics 
and quantum mechanically, see Refs. [38, 64, 68] and ref¬ 
erences therein. 

5 Bright exciton spin dynamics. Theory The dy¬ 
namics of bright exciton doublet can be conveniently 
described within the pseudospin formalism [64]. In this 
method, the 2x2 exciton spin-density matrix pK in the 
basis of 1^2 and Is states is decomposed as 

• a , (19) 

with nK = Tr{pK/2} being the spin-averaged distribu¬ 
tion function of excitons, S'k = Tr{pK^^/2} being the 
pseudospin. Here unit 2x2 matrix is omitted and a = 
{ax^ (Ty, a z) is 3. vector of Pauli matrices. Pseudospin com¬ 
ponents describe orientation of the oscillating microscopic 
dipole moment of exciton, particularly, Sz/n gives the de¬ 
gree of circular polarization or valley polarization, Sx/n 
and Sy/n give the degree of linear polarization in two axes 
frames rotated by 45° with respect to each other, i.e. valley 
coherences or exciton alignment. 

The exciton pseudospin satisfies the kinetic equa¬ 
tion [38,77] 

—^ X I7k = (20) 


27re^h‘^\(p{0)\^ {Kpcv){Kpl^) 

KijTriQUjQ K 


where I7 k is defined by Eq. (5) and QliSK} is the col¬ 
lision integral describing the variation of spin distribution 
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0 2 4 6 8 10 


Time, tjTc 

Figure 8 Examples of free exciton spin dynamics: Red 
solid curve shows the spin polarization as a function of 
time in the absence of scattering (r 2 ^ oo) calculated after 
Eq. (22) for Boltzmann excitons. Blue dashed curve shows 
exponential spin relaxation, Eq. (23), for finite scattering 
time r 2 /rc = 1/5; Tc is related with characteristic spin 
precession period in the field see text for details. 


function in the process of scattering. Its general form is 
given by the balance of in- and out-scattering processes: 

^ Sk' — H^k'k^'k) , (21) 

K' 


where the scattering rates IEkk' are determined by the par¬ 
ticular interactions, i.e. elastic scattering by impurities and 
other defects, exciton-phonon interaction, exciton-exciton 
and exciton-electron scattering. Equation (20) shows that 
the long-range exchange interaction induced effective field 
acts as a driving force for exciton spin dynamics. 

Similarly to the case of free excitons in quantum 
wells [38, 64] different regimes of spin decoherence in 
MX 2 monolayers can be identified depending on the rela¬ 
tion between the characteristic spin precession frequency 
and scattering rates [78]. To be specific, we address the 
situation where the excitons are assumed to be thermal- 
ized and distributed according to the Boltzmann function 
n{s) oc exjp {—s/ksT), where s = is exciton dis¬ 
persion, ksT is the exciton temperature measured in the 
energy units, ks is the Boltzmann constant. 

In order to illustrate the role of the effective field f2j^ 
in exciton spin/valley depolarization, we consider the hy¬ 
pothetic case where the scattering processes are completely 
neglected, = 0, but exciton distribution is Boltz¬ 

mann. This can be achieved if pump-induced initial dis¬ 
tribution is Boltzmann-like or if rapid exciton-exciton col¬ 
lisions quickly brought excitonic ensemble to some inter¬ 
nal thermal equilibrium and then become negligible due 
to the fast recombination process. Taking into account that 
the effective field is isotropic in the structure plane we ob¬ 
tain (neglecting the contribution of the small part of the 
wavevector space within the light cone) for thermalized ex¬ 


citons 


S,{t) = = Sz{0) 

K 


cos (l2K^)n(gK) 


^.(0) 




exp 



( 22 ) 


where erf ( 2 ;) = ^/2/n exp (—x^)dx is the error 

function, and — (^{kT)k‘^/h, with kT = 

^y2MkBT/W being thermal wavevector. According to 
Eq. (22) the exciton spin decay is accompanied with one 
half-period of oscillation, see red solid curve in Eig. 8. In 
this case the spin dephasing takes place only due to the 
spin precession: Excitons with different wave vectors have 
different precession frequencies resulting in the spin deco¬ 
herence. The decoherence rate I/tc can be estimated from 
the LT-splitting dX K ^ kT and scales as y/ksT refiecting 
the spanning of the LT splitting in the exciton ensemble. 
At long times, t Tc, one has strongly non-exponential 
decay, Sz{t)/Sz{0) ^ —r^/(2t^). The typical decay time 
of the in-plane pseudospin components Sx, Sy has the 
same order of magnitude as Tc. Equation (22) shows that 
even in the absence of scattering the long-range exchange 
interaction leads to the exciton spin/valley decoherence. 

Inclusion of scattering processes makes free exciton 
spin dynamics qualitatively different. In particular, in the 
strong scattering regime, where Out^sc ^ 1, where Tqc 
is the characteristic scattering time, the spin is lost by 
Dyakonov-Perel’ mechanism [79]. Hence, the spin decay 
law is exponential 


Si{t) = Si{Q)exp{-t/Tu), i = x,y,z, (23) 


see blue dashed curve in Eig. 8 and the components of spin 
relaxation rates tensor are given by [30, 38,64, 80] 

— = — = — = (24) 

'T~zz Tree 'T~yy 


where the angular brackets denote averaging over the en¬ 
ergy distribution and T 2 = T2(5k) is the relaxation time of 
second angular harmonics of the distribution function. This 
time describes also the decay of momentum alignment. In 
derivation of Eq. (24) elastic (or quasi-elastic) processes of 
exciton-defect and exciton-phonon scattering were taken 
into account in the collision integral, Eq. (21), in which 
case Wkk' = oc (5(£:k ~ ^k') ^^4 


1 

T2(eK) 


X^K'K(l-COs2t?'), 

K' 


(25) 


with being the angle between K' and x-axis. In the 
strong scattering regime, i 2 /c^r 2 ^ 1, the pseudospin ro¬ 
tation angles between the collisions are ^ 12 /c^T 2 ^ 1 and 
spin decoherence results from the random diffusion of the 
pseudospin over the sphere. It follows from Eq. (24) that 
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the temperature dependence of the relaxation rates is de¬ 
termined by an interplay of the increase of oc ksT 
with an increase in the temperature and possible energy, 
and thus temperature dependence of the scattering time T 2 . 

Equation (24) shows the anisotropy of exciton spin 
and valley relaxation. Particularly, circular polarization de¬ 
gree/valley polarization whose relaxation is governed by 
the rate 1/tzz decays twice faster compared with the lin¬ 
ear polarization degree/inter-valley coherence of excitons 
whose decay is controlled by smaller rate I/t^x = ^l^yy 
This is because the effective field due to the exchange inter¬ 
action lies in the structure plane, hence, both and 
f2y affect circular polarization degree, while linear polar¬ 
ization degree is affected by one component of only. 

Equation (20) is valid for free excitons, while experi¬ 
ments demonstrate that the localized exciton states play an 
important role in transition metal dichalcogenide monolay¬ 
ers [13, 18,23,40]. Eor completeness, let us discuss the spin 
decoherence of localized excitons, which can be caused by 
the fine structure splitting of localized electron-hole pair 
states provided that the localization potential has an in¬ 
plane anisotropy. The latter can be caused by a symmetry 
reduction induced by the coupling with the substrate, ran¬ 
dom deformations, anisotropic potential fluctuations, etc. 
We assume that the exciton is localized as a whole, i.e. its 
localization lengths a^' and a^/ along the main axes of x' 
and y' of the imperfection exceed by far excitonic effec¬ 
tive radius. Hence, in order to obtain an effective Hamilto¬ 
nian for a localized electron-hole pair one can average 1-Lx 
in Eq. (5) over the exciton center of mass wavefunction. 
As a result, we arrive at two eigenstates linearly polarized 
along the main axes x' and y' with the anisotropic splitting, 
cf. [68,77] 


AE, = 2 {{a{K)Kl,) - {a{K)K^,)) 


dyf Qjx 


(26) 


Here angular brackets denote quantum-mechanical aver¬ 
aging. The time scale for the spin decoherence of local¬ 
ized excitons is given by Tc, loc ^ kt/AEa. Very recently, 
the fine structure splittings for localized excitons were ob¬ 
served in WSe 2 monolayer flakes [81, 82]. The splittings 
are about 0.7 meV and are significantly larger than those in 
standard semiconductor quantum dots. This enhancement 
of the AEa is consistent with stronger excitonic effects, 
particularly, larger values of the Coulomb energy and of 
To. 

6 Bright exciton spin dynamics. Experiment The 

theory of bright exciton spin dynamics governed by the 
long-range exchange interaction outlined above has been 
developed for monolayers of transition metal dichalco- 
genides in Ref. [30]. It has been successfully applied to 
describe experimental data on MX 2 monolayers obtained 
by time-resolved photoluminescence on M 0 S 2 monolay- 




Figure 9 (a) Left axis: Time integrated photolumines¬ 
cence intensity as a function of emission energy. Right 
axis: Polarization of photoluminescence emission, (b) Left 
axis: photoluminescence emission intensity (black line) 
detected at maximum of A-exciton photoluminescence, 
E'det = 1.867 eV as a function of time and laser refer¬ 
ence pulse (green dashed line). Right axis: Calculated evo¬ 
lution of photoluminescence polarization after Eq. (23) us¬ 
ing Tzz ~ 4 ps (blue dotted line). T = 4 K, £’exc=1.965 eV, 

= 5 AtW//im2. Rgf [30] 


ers [30] and time-resolved Kerr rotation on WSe 2 mono- 
layers [23]. The role of the long-range exchange interaction 
in the intervalley excitonic transfer has been also confirmed 
by combination of time-resolved and cw polarized photolu¬ 
minescence spectroscopy of WSe 2 monolayers [24] and by 
cw photoluminescence spectroscopy on few-layer samples 
of M 0 S 2 [83]. 

Eigure 9 presents the results of photoluminescence ex¬ 
periments carried out in Ref. [30] on M 0 S 2 monolayer 
deposited on the Si02/Si substrate. These experiments as 
well as the experiments in Ref. [23] and described be¬ 
low have been performed with the samples in vacuum. 
Both photoluminescence intensity and its circular polar¬ 
ization degree (under circularly polarized excitation) were 
recorded as a function of emission energy. As a simplest 
possible model, we assume that the stationary, i.e. time in¬ 
tegrated, polarization is determined by the initially created 
polarization Pq, the lifetime of the electron-hole pair r and 
the polarization decay time as [84]: 


P.= 


1 + 't/'Ts 


(27) 


In Eig. 9(a) an average, time-integrated photoluminescence 
polarization of Pc ~ 60% in the emission energy range 
1.82 ... 1.91 eV is observed. The emission time measured 
in time-resolved photoluminescence experiments is r 
4.5 ps, extracted from Eig. 9(b) in agreement with earlier 
measurements [20]. Eor the studied M 0 S 2 it was not pos¬ 
sible to conclude if the measured emission time is a radia¬ 
tive lifetime or it is limited by non-radiative processes. Eor 
Po = 100% we find an estimate of 7 ps. This value is 
in reasonable agreement with crude theoretical estimates 
of Tzz^^ ps made for this sample under the assumption 
that the spread of excitons in the energy space is limited by 
the collisional broadening, ^ /i/r 2 , rather than by the ki- 
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Figure 10 Kerr rotation dynamics at T = 4 K for cr+ 
and CF~ pump beam. The laser excitation energy is Ei = 
1.735 eV. Inset shows sketch of time resolved Kerr rota¬ 
tion effect. Data from Ref. [23]. 


netic energy distribution [30,38]. Note, that in this case the 
criteria of applicability of kinetic Eq. (20) are not fulfilled, 
which results in inaccuracy of Tzz estimation. It is worth 
stressing, that even with a value of Tzz in the picosecond 
range a high polarization can be obtained in time integrated 
measurements. This is because the photoluminescence de¬ 
cay time r is also ultrashort. It is demonstrated in Fig. 9(b), 
which allows to compare the theoretically predicted polar¬ 
ization decay and the measured photoluminescence inten¬ 
sity decay. 

Time-resolved Kerr rotation technique serves as an al¬ 
ternative and highly sensitive tool to study spin dynamics, 
see Refs. [86,87] for review. In this method the spin-Kerr 
effect, i.e. polarization plane rotation angle of the linearly 
polarized probe pulse, 0{At), is measured as a function of 
the time delay. At, between the probe and preceding circu¬ 
larly polarized pump pulse, see inset in Fig. 10. The pump 
pulse creates spin/valley polarization of excitons, which re¬ 
sults in polarization-dependent modification of the param¬ 
eters cjQ, Fq and F in the MX 2 reflection coefficient, r(u;), 
Eq. (8) due to nonlinearities in the system [87]. Thus, the 
probe beam monitors imbalance of cr+ and a~ excitons in 
the system. 

Such measurements were carried out in Ref. [23] on 
high quality WSe 2 monolayers where the trion and exciton 
lines in photoluminescence are separated by ~ 30 meV and 
the exciton line is stable up to the room temperature. The 
latter makes it possible to experimentally address tempera¬ 
ture dependence of exciton spin dynamics, unaccessible in 
conventional semiconductors. Additional advantage of the 
Kerr rotation technique comes from the fact that the limited 
time resolution of photoluminescence spectroscopy does 
not allow one to address spin dynamics at elevated temper¬ 





Figure 11 (a) Kerr rotation dynamics after a cr+ polarized 
pump pulse for different lattice temperatures; (b) Transient 
reflectivity dynamics for different temperatures. The laser 
excitation energies are identical to the ones used in (a), see 
text; (c) Temperature dependence of the measured (sym¬ 
bol) and calculated after Eq. (24) (solid and dashed lines) 
exciton valley polarization relaxation time, see text for de¬ 
tails. From Ref. [23]. 


atures where spin relaxation time decreases [20,40] . Kerr 
rotation dynamics has higher time resolution (^100 fs) and 
makes it possible to apply a strictly resonant excitation of 
the exciton. 

Figure 10 shows the Kerr rotation dynamics, 0{At), 
measured at 4 K for both cr+ and polarized pump 
pulses. The pump energy Ei = 1.735 eV is set to the max¬ 
imum of the Kerr signal, which is very close to the neu¬ 
tral exciton transition identified in the photoluminescence 
spectra. The observed sign reversal of the Kerr signal in 
Fig. 10 at the reversal of pump helicity is a consequence 
of the selective optical initialization of the and K_ 
valley, respectively. The transient reflectivity measured at 
the same conditions using linearly cross-polarized pump 
and probe pulses is shown in Fig. 11(b). It is seen that the 
reflectivity decay time is about ten times longer than the 
one observed in time-resolved Kerr rotation. Therefore the 
mono-exponential decay time = (6 dz 0.1) ps of the 
Kerr rotation dynamics at T = 4 K in Fig. 10, probes di¬ 
rectly the fast exciton valley depolarization. The value of 
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Ts is close to that measured by time-resolved photolumi¬ 
nescence on M0S2 samples, see above. 

Figure 11(a) displays the variation of the spin-Kerr ef¬ 
fect dynamics as a function of the temperature. For T > 
30 K, a clear decrease of the exciton spin/valley polar¬ 
ization decay time is observed down to 1.5 ps at T = 
125 K. At the same conditions reflection coefficient decay 
is much longer, see Fig. 11(b). The excitation power de¬ 
pendences of the exciton dynamics have been also investi¬ 
gated. In the studied power range corresponding to varia¬ 
tion of exciton density from ~ 1.5x 10^^ to 10^^ cm 
both the Kerr rotation and reflectivity dynamics do not de¬ 
pend on the photo-generated exciton density within the ex¬ 
perimental accuracy. This demonstrates that the exciton- 
exciton interactions play a minor role in exciton valley dy¬ 
namics presented in Figs. 10 and 11 over this density range. 
Note, that in Ref. [24] much higher excitation intensities 
were applied to WSe 2 monolayers resulting in the exci¬ 
ton densities up to Ux ^ 10^^ cm“^. In this case the ef¬ 
fect of exciton density on circular polarization of emission 
was observed. Particularly, Ref. [24] reports an increase of 
the circular polarization degree with an increase of exciton 
density up to 5 x 10^^ cm“^ and a decrease of the polar¬ 
ization degree with further increase in Ux. 

Figure 11(c) shows the experimentally measured in 
Ref. [23] exciton polarization decay times and theoretical 
calculations carried out after Eqs. (14) and (24) for n = 
\/T0 and Jo = 0.16 ps“^. Qualitatively, the drop of the 
exciton spin/valley relaxation time when the temperature 
increases can be well explained by the increase with the 
temperature of the precession frequency I7 k oc K, which 
makes spin precession and decoherence faster, Eq. (24). 
The solid and dashed lines in Fig. 11(c) correspond to the 
calculated exciton spin/valley relaxation time for two dif¬ 
ferent scattering times T 2 = 0.066 ps and T 2 = 0.25 ps, 
respectively. 

The choice of the scattering time values in the calcu¬ 
lations of Ref. [23] was motivated by the following rea¬ 
sons. The latter (longer) value of T 2 corresponds to exci¬ 
ton energy uncertainty equivalent to 30 K corresponding 
to loffe-Regel criterion of delocalization where the prod¬ 
uct of the thermal wavevector kr and the mean free path I 
is on the order of 1. However, a much better agreement 
between the calculated and measured exciton relaxation 
times in Fig. 1 1(c) for T > 30 K is observed for the smaller 
value of the scattering time T 2 = 0.066 ps, but formal cri¬ 
terion of kinetic equation is fulfilled in this case only for 
high temperatures, T > 100 K. The possible origins of the 
scattering are: (i) short-range defects, (ii) Coulomb scat¬ 
tering with resident electrons. In the latter case for resident 
electron density ^ 10^^ cm“^ the exchange exciton- 
electron [88] scattering yields T 2 ^ 10“^^ s. 

At low temperatures, 4 < T < 30 K, the measured ex¬ 
citon spin relaxation time is almost temperature indepen¬ 
dent. This behavior could be either caused by (i) the very 
strong scattering in a regime where is smaller than the 
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Figure 12 (a) TRKR traces measured on the M 0 S 2 flake 
for different helicities of the pump beam, with a laser ex¬ 
citation energy of 1.937 eV. (b) Normalized time resolved 
Kerr rotation (TRKR) and transient reflectivity traces mea¬ 
sured on the same flake with a laser excitation energy of 
1.937 eV. The solid line indicates a biexponential fit to the 
time resolved Kerr rotation trace. The inset shows a high- 
resolution plot of the data. The solid line indicates a mo¬ 
noexponential fit to the transient reflectivity trace. From 
Ref. [89]. 


collision broadening leading to a temperature independent 
spin relaxation time, see Refs. [30, 38] and discussion of 
experiment on M 0 S 2 above, or (ii) by a localized character 
of the exciton below 30 K. Valley dynamics of excitons in 
WSe 2 has been also addressed by photoluminescence spec¬ 
troscopy [24] and by Kerr rotation spectroscopy [85]. The 
origin of discrepancy of the results of Refs. [24] and [85] 
as well as of Ref. [23] is not clear and is probably related 
to different experimental methods. 

In experiments performed in Ref. [23] no long-living 
spin/valley polarization was detected: The Kerr rotation 
signal is absent for the pump-probe delays At ^ 25 ps, 
where photoinduced reflectivity vanishes. It indicates ab¬ 
sence of polarization transfer from excitons to resident 
carriers in WSe 2 , unlike similar time-resolved Kerr ex¬ 
periments performed in III-V or II-VI semiconductor 
nanosystems, see Refs. [86, 87] for review. By contrast, 
experiments performed on M 0 S 2 monolayer flakes in 
Ref. [89] indicate valley polarization of resident carri¬ 
ers. Figure 12(a) demonstrates Kerr rotation signal as a 
function of the pump-probe delay for and ( 7 ~ helic¬ 
ities of the pump pulse. Like the data shown in Fig. 10 
the signal flips its sign with reversal of radiation helicity 
evidencing selective population oi K± valleys. However, 
unlike data shown in Fig. 11, the transient reflectivity mea¬ 
sured in Ref. [89] on M 0 S 2 flake decays on the timescale 
of ^ 4.5 ps, while the Kerr rotation signal demonstrates 


Copyright line will be provided by the publisher 















pss header will be provided by the publisher 


13 


bi-exponential behavior, see Fig. 12(b). The analysis car¬ 
ried out in Ref. [89] shows that this decay contains two 
components, the fast one with the same decay time as that 
of the reflectivity, and the slow one with the longer de¬ 
cay time of 17 ps at T = 4.5 K. This long decay time 
is attributed in Ref. [89] to the spin/valley polarization of 
resident electrons in the sample. In this case time resolved 
Kerr rotation technique provides an access to the resident 
electron spin dynamics in transition metal dichalcogenide 
sample. 

7 Outlook The understanding of spin and valley dy¬ 
namics of excitons in transition metal dichalcogenides has 
rapidly evolved over the last three years. Although this area 
of research started only in 2012, a number of solid ex¬ 
perimental facts has been already established and conclu¬ 
sive theoretical models have been developed. Particularly, 
it has been demonstrated that the spin and valley dynamics 
of bright exciton doublet in MX 2 monolayers is governed 
by the long-range exchange interaction between an elec¬ 
tron and a hole forming an exciton. This interaction acts, 
in the pseudospin picture, as an effective magnetic held 
coupling the circularly polarized states of moving exciton. 
The long-range exchange interaction can also be responsi¬ 
ble for the localized excitons flne structure and their spin 
decoherence. 

Still, further progress in needed before a complete un¬ 
derstanding of multifaceted excitonic spin dynamics in 
MX 2 monolayers can be reached. Despite strong theoret¬ 
ical and experimental progress the full picture of the spin 
dynamics of charged excitons is still absent. It is still not 
clear why no optical orientation of excitons is observed 
in high-quality MoSe 2 monolayers [90]. Moreover, the 
magnetooptical studies of excitons and trions in transition 
metal dichalcogenides such as WSe 2 and MoSe 2 reveal 
large exciton ^-factors in magnetic held perpendicular to 
the sample plane, \gx\ ^ 4 [52,91-94]. These observa¬ 
tions are very surprising because the two-band k • p model, 
Eq. (1) used here, results in = 0 [52]. Therefore fur¬ 
ther experimental and theoretical studies are needed to 
elaborate the band structure description of transition metal 
dichalcogenides monolayers. 
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